Ground-state phase diagram of the one-dimensional Hubbard model with an 

alternating potential 
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We investigate the ground-state phase diagram of the one-dimensional half-filled Hubbard model 
with an alternating potential — a model for the charge-transfer organic materials and the ferro- 
electric perovskites. We numerically determine the global phase diagram of this model using the 
level-crossing and the phenomenological renormalization-group methods based on the exact diago- 
nalization calculations. Our results support the mechanism of the double phase transitions between 
Mott and a band insulators pointed out by Fabrizio, Gogolin, and Nersesyan [Phys. Rev. Lett. 83, 
2014 (1999)]: We confirm the existence of the spontaneously dimerized phase as an intermediate 
state. Further we provide numerical evidences to check the criticalities on the phase boundaries. Es- 
pecially, we perform the finite-size-scaling analysis of the excitation gap to show the two-dimensional 
Ising transition in the charge part. On the other hand, we confirm that the dimerized phase survives 
in the strong-coupling limit, which is one of the resultants of competition between the ionicity and 
correlation effects. 

PACS numbers: 71.10.Pm, 71.30.+h 



I. INTRODUCTION 

The electronic and/or magnetic properties of the low- 
dimensional interacting electrons have attracted great 
interest in researches of materials, such as the quasi 
one-dimensional (ID) organic compounds and the two- 
dimensional (2D) high-T c cuprates, where a variety of 
generalized Hubbard-type models have been introduced^ 
For the ID case, a concept of the Tomonaga-Luttingcr 
liquid (TLL) has been widely accepted and intensively 
used not only for the descriptions on the low-energy 
and long-distance behaviors of the critical systemsj^*^ 
but also for the prediction of its instabilities to, for in- 
stance, various types of density-wave phases observed in 
the models^ 

The ID Hubbard model with an alternating potential 
(also called the ionic Hubbard model) is one of the mod- 
els for the 7r-electron charge-transfer organic materials, 
such as TTF-Chloranil|& and/or the ferroelectric tran- 
sition metal oxides as BaTiOsA^ It is defined by the 
Hamiltonian 



C j,s C j+l,s + H.c. 
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5>(-l)^, (1) 



where Cj a annihilates an s-spin electron (s =| or !) on 
the jth site and the number operator n,j tS — cj s Cj s and 
Uj = Tijrf + nj,i- While t and U terms stand for the elec- 
tron transfer among sites and the Coulomb repulsion on 
the same site, respectively, the A term represents an en- 
ergy difference between the donor and acceptor molecules 
(or between the cation and oxygen atoms), and it intro- 
duces ionicity effects into the correlated electron systems 
(wc set t = 1 in the following discussion). 



The understandings on the model have been accumu- 
lated in the literature, where the theoretical investiga- 
tions including numerical calculations have been per- 
formed mainly at the half filling: Nagaosa and Taki- 
moto calculated the magnetic and charge-transfer gaps 
as functions of A (U fixed) by using the quantum Monte 
Carlo (QMC) simulation^ Resta and Sorella, using the 
exact-diagonalization calculations of finite size systems, 
reported, for instance, the divergence of the average dy- 
namical charged By applying the renormalization-group 
(RG) method to the bosonized Hamiltonian, Tsuchiizu 
and Suzumura estimated a boundary line between the 
Mott insulator (MI) and a band insulator (BI) phases in 
the weak-coupling regions^ On the other hand, Fabrizio, 
Gogolin, and Nersesyan (FGN) predicted an existence 
of the "spontaneously dimerized insulator" (SDI) phase 
between them^SiAi After their proposal, various numeri- 
cal calculation methods have been so far applied to con- 
firm it: Wilkens and Martin performed the QMC simu- 
lations to evaluate, e.g., the bond order parameter, and 
reported the transition between the BI and SDI phases 
and stated an absence of MI phase for A > Oii^ By the 
combined use of the method of topological transitions 
(jumps in charge and spin Berry phases and the 
method of crossing excitation levels, Torio et al. pro- 
vided a global ground-state phase diagram, which is in 
accord with the FGN scenario^ And an existence of 
the SDI phase for all U > regions was first exhib- 
ited there. The density matrix renormalization group 
(DMRG) calculationsiiiSiiiLSi have been performed by 
several groups. For instance, Zhang et al. provided the 
data on the structure factors of relevant order parame- 
ters in the weak- and intermediate-coupling region, which 
supports an existence of intermediate SDI phase between 
the BI and MI phases^ On one hand, Kampf et al. es- 
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timated the excitation gaps up to 512-site system and 
found the boundary of the BI phase while the existence 
of the second boundary was not resolvediSi Therefore, 
some controversy as well as points of agreement exists in 
these recent investigations. 

In this paper using the standard numerical techniques, 
we shall provide both the global structure of the ground- 
state phase diagram and the evidences to show the crit- 
icalities of the massless spin and charge parts. For this 
purpose, it is worthy of noting that the FGN scenario 
consists of two types of instabilities commonly observed 
in the TLL, i.e., the transition with the SU(2)-symmetric 
Gaussian criticality in the spin part, and that with the 
2D-Ising criticality in the charge part (see Sec.[n]). Fur- 
thermore, these types of phase transitions have been 
numerically treated by the level-crossing (LC) method, 
and the phenomenological renormalization-group (PRG) 
methodiSi The LC method has been applied to the frus- 
trated XXZ chain;22i2i and also used in the research 
of higher-S 1 spin chains^ spin ladders^ and ID cor- 
related electron systemsiSSiSi The advantage of using the 
LC method is not restricted to its accuracy in estimat- 
ing the continuous phase transition points including the 
Berezinskii-Kosterlitz-Thouless type one; it also provides 
a means to check their criticalities (see Sec. IIII|) . 23 Both 
of these are important in order to settle the controversy 
mentioned above, and, in fact, the precise estimation of 
the spin-gap transition point of the S = | J1-J2 chain 
was first given by the LC method^ while numerical in- 
vestigations including the DMRG work were performed. 
On the other hand, the PRG method is also a reliable nu- 
merical approach to determine second-order phase tran- 
sition point, especially for the 2D-Ising transition where 
the LC method is not available. Analysis based on the 
PRG method for the 2D-Ising transition is successful 
in the spin systems^* Furthermore, one of the authors 
treated the 2D-Ising transition in the S = | J1-J2 model 
under a staggered magnetic field, where the critical phe- 
nomena in the vicinity of the phase boundary line were 
argued^ Therefore, based on these recent developments, 
we shall perform the numerical calculations; especially, to 
our knowledge, this is the first time that the PRG method 
successfully applied to the 2D-Ising transition observed 
in one part of the two-component systems like the inter- 
acting electrons. 

The organization of this paper is as follows. In Sec.lTTl 
we shall briefly refer to the effective theory based on the 
bosonized Hamiltonian and order parameters of expected 
density-wave phases, and mention the FGN scenario. In 
Sec. IIIII we explain procedures of the numerical calcula- 
tion to determine transition lines, where connections be- 
tween the methods and instabilities of the TLL systems 
will be explained. After that, we shall give a ground-state 
phase diagram in whole parameter region. Furthermore, 
to confirm the criticalities and to serve a reliability of our 
calculations, we check the consistency of excitation lev- 
els in finite-size systems. A finite-size scaling analysis of 
the charge excitation gaps is also performed in the vicin- 



ity of the phase boundary line. Section IIVI is devoted 
to discussions and summary of the present investigation. 
A short comment on the Berry phase method2iiiiiiiSii& 
will also be given there. We will provide the comparison 
with that method, which is helpful to exhibit a reliability 
of our approach as well as the results. 



II. GROUND STATES AND PHASE 
TRANSITIONS 

The bosonization method provides an efficient way 
to describe low-energy properties of the ID quantum 
systems^ Linearizing the cos-band at two Fermi points 
±fcp = ±7rn/2a [an electron density n := N/L = 1 and a 
number of sites (electrons) L (N)}, and according to stan- 

is given 



dard procedure, the effective Hamiltoniar-^ 
as H -> H = H„ + H a + H 2 with 



H v 



H 2 



dx ~r- 
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yfajh, (v = p,a), (2) 
sin V2<j) P cosV2> CT . (3) 



The operator V is the dual field of <j> v satisfying the com- 
mutation relation [4> u (x), d y 9 u i{y) /-k] = iS(x — y)8 v , v > ■ K v 
and v v are the Gaussian coupling and the velocity of el- 
ementary excitations. Coupling constants g p (< 0) and 
g a stand for the 4fcp-Umklapp scattering and the back- 
ward scattering bare amplitudes, respectively, and Hi ex- 
presses a coupling between the spin and charge degrees 
of freedom. In Table we summarize the order param- 
eters for the relevant 2fcp density- wave phases, i.e., the 
charge-density- wave (CDW), bond charge-density- wave 
(BCDW), and spin-density-wave (SDW) phases, where 
the electron's spin and the bond charge are given as Sj = 

E s ,y 4, s[5 CT Ws>' and n i = E«( c J,» c i+i.« + H - c -)' re ~ 

spectively (cr are the Pauli matrices). Their bosonized 
expressions are given in the second column. In the third 
column, we give the locking points of phase fields. As 
discussed in Ref. fToL there arc two locking points of <f> p , 
i.e., (y/8(f> p ) = ±0o in the BCDW state. The phase 4>o, a 
function of U and A, continuously varies from to n. 

Let us see the system with increasing A for fixed U. 
At A = 0, the ground state is in the MI phase with 
the most divergent SDW fluctuation (the third row of 
Table HJ). According to the arguments the MI phase 
may survive for U 2A. For 2A U, H2 becomes 
relevant, and leads to the BI phase with the long-range 
CDW order without degeneracy (the first row). For this 
issue, FGN argued that under the uniform charge distri- 
bution a renormalization effect of H2 to g a brings about 
the spin-gap transition in the spin part at a certain value 
of A a (U), which is described by the sine-Gordon (SG) 
theory. This is qualitatively in accord with the perturba- 
tion calculation in the strong-coupling region^ and leads 
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TABLE I: The order parameters. The bosonized forms and 
the locking points of phase variables ({\/8(/) p }, (\f&<f> a )) are 
given in the second and third columns. (f>o is a function of U 
and A, and * denotes a phase not to be locked. 

Order parameters Bosonized forms Locking points 
Ccdw = ( — 2 sin V2(f> p cos y/2cj> a (tt,0) 
Obcdw = (— 2 cos \/20p cos \/20<t (±0 O ,O) 



to the SDI phase with the long-range BCDW order (the 
second row). Further with the increase of A, a tran- 
sition in the charge part occurs on a separatrix A P (U) 
between two different types of charge-gap states. This 
line corresponds to the massless RG flow connecting the 
Gaussian (the central charge c = 1) and the 2D-Ising 
(c = |) fixed points^ and its description is given by the 
double-frequency sine-Gordon (DSG) theory^ Our main 
task is thus to estimate A U (U) for U > and to check 
the criticalitics based on their prediction. 

III. NUMERICAL METHODS AND 
CALCULATION RESULTS 

Low-lying excitations observed in the finite-size sys- 
tems are expected to serve for the determinations of tran- 
sition points. Here, we take a look at the following oper- 
ators with lower scaling dimensions: 

0„,i = V2cos%/2c/v, (4) 
0„,a = \/2sin\/2<^, (5) 
0„ )3 = exp(±iV20„). (6) 

According to the finitc-size-scaling argument based on 
the conformal field theory, corresponding energy levels 
for these operators AE V ^ (taking the ground-state en- 
ergy Eq as zero) are expressed by the use of their scaling 
dimensions i^jifi 

2ttv u 

AE v . l ~ — —x Vti . (7) 

Then these excitations can be extracted under the 
antipcriodic boundary condition with respect to the 
ground state due to the selection rule of the quan- 
tum numbers iSS*21 In the numerical calculations using 
the Lanczos algorithm we can identify AE V ^ according 
to the discrete symmetries of the wave functions, e.g., 
translation (cj <a — > Cj+2, s ), charge conjugation [cj >3 — > 
(— l) J 'ct +1 ], spin reverse (cj jS — > Cj t _ a ), and space in- 
version (cj, s — » CL-j, s )- Here note that, except for the 
spin-reversal operation, definitions of these transforma- 
tions arc different from those of the uniform systems, 
such as the extended Hubbard model M- 

First, we tre at the spin -gap transition in the spin part 
following Refs. I22l2& and[27[ In the SDW phase, due to 



the marginal coupling in the SU(2)-symmetric spin part, 
the singlet (av,i) and triplet {x a ^ = x a ^) excitations 
split as x a} i > x a} 2 = £<r,3 satisfying a universal relation 

Xa-i + 3x ai 2 1 fo . 

4— 2" (8) 

Then, the degeneracy condition 

x a ,i = x a<2 = x ai3 (9) 

stands for the vanishing of the coupling, and provides a 
good estimation of the spin-gap transition point. Note 
that Torio et al. used the crossing of these excitation 
levels for the determination of the MI-SDI transitions- 
while the consistency check of the levels to confirm 
the universality of transition is still absent. Figure ^ 
shows an example of the A dependences of x a .i for the 
16-site system at u = 0.6 [here we introduce the re- 
duced Coulomb interaction parameter u = U/(U + 4)]. 
For this plot, we estimated the spin- wave velocity v a 
from a triplet excitation with the wave number An/L as 
v a = linii^oo AE(S = l,k = 4m / L) / {2-k / L) and nor- 
malized the excitation gaps AE a> i according to Eq. (7J). 
The singlet (triplet) level corresponding to the operator 
Co-,1 [0(7,2 (0(7,3)] is denoted by circles (triangles) with 
a fitting curve. Their behaviors reflect the TLL prop- 
erties: For instance, the amplitude of the level splitting 
decreases with the increase of A due to its renormaliza- 
tion effect, and eventually the level crossing occurs at 
A cr (/7, L). More precisely, in order to confirm the univer- 
sality, we plot the averaged scaling dimension x av , i.e., 
the left-hand side of Eq. (JSJ in Fig. ^ (squares). We also 
exhibit the L dependence of x av at A = 1 .0 as an exam- 
ple (see the inset). The result shows that the condition 
imposed on x a ^ is accurately satisfied for A < A a (U, L); 
in particular, the extrapolated value of x av is almost i. 
Consequently, the level crossing at which Eq. © is sat- 
isfied can be regarded as an indication of the spin-gap 
transition in the spin part of the Hamiltonian . On the 
other hand, the spin part is dimerized for A > A a (U, L). 

Next, we discuss the 2D-Ising transition in the charge 
part. Recently, we have treated the crossover behavior 
into the 2D-Ising criticality in the study of the frustrated 
quantum spin chain^ so we shall here employ the same 
approach to determine A P (U). Since there are two criti- 
cal fixed points connected by the RG flow, a relationship 
between lower-energy excitations on these fixed points 
is quite important. For this, the so-called ultraviolet- 
infrared (UV-IR) operator correspondence provides sig- 
nificant informations iiLiS Along the RG flow, the oper- 
ators on the Gaussian fixed point (UV) are transmuted 
to those on the 2D-Ising fixed point (IR) as 

O pA ^n, p , 2 ^/ + e, (10) 

where fi is the disorder field (Z2 odd), and e is the en- 
ergy density operator (Z2 even) with scaling dimensions 
x p = g and x e = 1, respectively. Furthermore, since a 
deviation from the transition point A — A p (U), which 
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FIG. 1: The A dependence of x a ,i at u = 0.6 for the 16- 
site system [u — U/(U + 4)]. The spin-gap transition point 
Ao-(t/, L) is estimated from the level crossing between the 
singlet (circles) and triplet (triangles) spin excitations. The 
squares plot x av = {x a ,i + 32V,2)/4, and the inset shows the 
L dependence of a; av at A = 1.0, where a least-square-fitting 
line to the data of L = 12-16 is given. 



is the coupling constant of the p> 2 term in the DSG 
Hamiltonian^i plays a role of the thermal scaling vari- 
able, anomalous behaviors in the vicinity of A p (U) are to 
be related to the divergent correlation length of the form 
£ oc [A - A P {U)]~ V with the exponent 1/v = 2 - x e = 1. 
On one hand, the excitation [i corresponding to Pi i pro- 
vides a lower-energy level, so we shall focus our attention 
on it. In order to determine the transition point, we shall 
numerically solve the following PRG equation for a given 
value of U with respect to Ai2L2& 

(L + 2)AE Ptl {U,A,L + 2)=LAE Ptl (U,A,L). (11) 

Since this is satisfied by the gap AE P: i (U, A, L) oc 1/L, 
the obtained value can be regarded as the L-dependent 
transition point, say A p (U,L + 1). We plot L and A 
dependences of the scaled gap LAE Pt i(U, A, L) in Fig. EJ 
and find that the size dependence of the crossing point is 
small for large values of U, but it is visible in the weak 
coupling case. 

While the results in the thermodynamic limit will be 
given in the last part of this section, we shall check first 
the criticality on and in the vicinity of the phase bound- 
ary using the extrapolated data A p (U). For this aim, an 
evaluation of the central charge c through the size depen- 
dence of the ground-state energy provides a straightfor- 
ward way^ However, as exhibited in the following, the 
critical line in the charge part is close to the spin-gap 
transition line, so that influences from the spin part with 
the small dimer gap prohibit a reliable estimation of c 
from the data of the finite-size systems. Alternatively, 
we shall evaluate a ratio of the charge-excitation gaps 
AE p>1 (U, A, L) and AE Pt2 (U, A, L) on the phase bound- 




FIG. 2: The L and A dependences of the scaled gap LAE p a. 
From left to right, u =0.12, 0.60 and 0.72, respectively. The 
correspondence between marks and system sizes is given in 
the figure. Crossing points give the L-dependent transition 
points A P (U,L + 1). 

ary to check the UV-IR operator correspondence. Ac- 
cording to Eqs. (J7J) and lfTU|l . it is expressed by the scaling 
dimensions of operators e and fj, as 

AE pA (U,A p (U),L) x E= l 
AE Pt2 (U,A p (U),L) ~" x e 8 [ ' 

for large L. Figure [3] plots the A dependence of R for 
L = 10-16 (u = 0.72). The transition point in the ther- 
modynamic limit is denoted by the arrow near the x axis. 
While the ratio exhibits a subtle A dependence around 
the point, we interpolate these data, and estimate the L 
dependence of R at A p (U), which is given with a least- 
square-fitting line in the inset. The plot shows that the 
extrapolated value is fairly close to | . Therefore we con- 
clude that the boundary line A p (U) belongs to the 2D- 
Ising universality class. 

Furthermore, we shall investigate the critical 
behavior*^ According to the finite-size-scaling ar- 
gument, we analyze the charge-excitation gap by using 
the following one-parameter scaling form: 

AE pA (U, A, L) = L-^(L[A - \(U)] V ). (13) 

Since AE p i cx l/£ in the thermodynamic limit (L/£ — > 
oo), the scaling function is expected to asymptotically 
behave as ty(x) oc x for large x. On the other hand, the 
gap AE Pt i oc 1/L on the critical point (£/£ — * 0) so that 
ty(x) ~ const for x — > 0^ Figure 0] plots Eq. (fTT^ft using 
the exponent of the 2D-Ising model v = 1. Although due 
to the smallness of L a scattering of the scaled data is vis- 
ible especially near the transition point, the data of dif- 
ferent system sizes are collapsed on the single curve, and 
its asymptotic behaviors agree with the expected ones. 
Therefore, we can check that, in the transition of the 
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FIG. 3: The A dependence of the charge-excitation-gap ratio 
R = AE p , 1 (U,A,L)/AE Pz2 (U,A,L) for L = 10-16 at u = 
0.72. The arrow shows the transition point A P (U). The inset 
plots the L dependence of R at A P (U) with a least-square- 
fitting line. 
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so here we perform the same calculations in order to 
complete the ground-state phase diagram. We employ 
the formula: A a (U,L) = A a (U) + aL~ 2 + bL~ A , where 
A a (U), a and b are determined according to the least- 
square-fitting condition. Then, we extrapolated the data 
of L = 12-18 as shown in Fig. 0Ja), where from bot- 
tom to top the data with fitting curves are given in the 
increasing order of U. Consequently, the spin-gap tran- 
sition line A cr (C/) (open circles with a fitting curve) is 
given in Fig. where the reduced alternating potential 
parameter S = A/(A + 2) is used as the y axis. On 
the other hand, for the extrapolation of A p (U, L), we as- 
sume the following formula^ A P (U, L) = A p (U) + aL~ 3 , 
and extrapolate the data of L = 10-18 as shown in 
Fig. [Hfb). Consequently, Fig. [3] shows that the criti- 
cal line in the charge part (open squares with a fitting 
curve) does not coincide with the spin-gap transition line, 
i.e., A a (U) < A P (U), and that the 2D parameter space 
{{u,6) | < u, 5 < 1} is separated into the MI, BI, and 
SDI phases with SDW, CDW, and BCDW, respectively. 
Since the Hubbard gap provides a principal energy scale 
and a shape of the boundary is roughly determined so 
that the magnitude of the band gap becomes compara- 
ble to the scale, the U dependence of the boundaries is 
expected to be weak in the small-?/ regionj&iS which is 
in agreement with our observation. On the other hand, 
in order to clarify the behaviors in the large-?/ region, 
we plot a magnification of the phase diagram around the 
2A = U line in Fig. |SJ This shows that in the limit of 
U — > oo the boundaries do not merge to the line: More 
precisely, for U — 96 we obtain A p — U/2 ~ —0.65 and 
A a - U/2 ~ -0.97, respectively. In Rcf. H adding to 
the spin part (2A,j -U ~ -1.91 for U, V > 1), they also 



reported 2A^ 



U 



1.33, which is close to our esti- 



mation. Consequently, we confirm that the intermediate 
SDI phase may survive in the large-C/ limit, which is one 
of the nontrivial behaviors and is contrasted to the naive 
argument. 



FIG. 4: The finite-size-scaling plots of the charge-excitation 
gap AE p ,i for systems of L = 14-18 at u = 0.72 and 0.80. We 
use the 2D-Ising critical exponent v — 1. A dotted line (the 
slope 1) is given for the guide to eye. 



charge part, the deviation A — A p (U) plays a role of 
the thermal scaling variable on the 2D-Ising fixed point. 
Here, note that in the strong-coupling region the energy 
scale of the crossover behavior may be large enough to 
be detected even in the small-size systems. However, the 
finite-size-scaling nature may become obscure in the weak 
and intermediate couplings. 

Lastly, we present the ground-state phase diagram. In 
order to determine it, the extrapolations of A V (U,L) to 
the thermodynamic limit are carried out. For the spin 
part, it should be noted that Torio et al. evaluated the 
spin-gap transition line from the level crossing Eq. JJJJ, 16 



Here we shall perform a comparison with the previ- 
ous DMRG results. As mentioned in Sec. while the 
DMRG calculations performed by several groups seem 
not to reach an agreement with respect to an existence 
of the SDI phase, it may be informative to provide a 
comparison with our result. Zhang et al. determined 
two-types of phase transition points U c \ and U C 2 based 
on the structure factor of the BCDW order parameter^ 
we plot their results in Fig. O by using the filled squares 
and filled circles, respectively. This shows that their esti- 
mations of U c i agree well with our data A P (U), although 
those of U C 2 considerably deviate from A a (U). Since 
the phase transition at A a {U) is the spin-gap transition, 
the logarithmic corrections to the power-law behaviors 
as well as the exponentially small magnitude of the spin 
gap generally make it difficult to determine the transi- 
tion point. On one hand, as explained in the above, the 
LC method used here overcomes these difficulties in the 
determination of the transition points A CT (/7). 
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FIG. 5: The ground-state phase diagram of the ID Hub- 
bard model with the alternating potential. The open circles 
(squares) with a fitting curve show the spin-gap (2D-Ising) 
transition line in the spin (charge) part. The stable regions 
of the MI, SDI, and BI phases are given in the 2D parameter 
space (u, 5) [u = U/(U + 4) and 8 = A/(A + 2)]. Insets (a) 
and (b) show the extrapolations of the //-dependent transi- 
tion points in the spin and the charge parts, respectively. For 
comparison, we also plot the DMRG calculation results given 
in Ref. by using the filled squares (U c i in their notation) 
and the filled circles {U C 2)- 



IV. DISCUSSION AND SUMMARY 

For the understanding of the phase diagram in the 
largc-[7 limit, let us see the perturbativc treatment of 
Hamiltonian Q under the condition of U — 2 A ^> 1. 
A a {U) may be related to the spin-gap transition point in 
the S = i J\- J2 modeliA2*2S Therefore, using its numer- 
ical valued and perturbative expressions on J\ and Ji^ 
we can approximately estimate A a (U) as a solution of 
the equation J 2 /J 1 ~ X/(l - AX) ~ 0.2411, where X = 
(l + 4x 2 -x 4 )/U 2 (l-x 2 ) 2 andx = 2A/U. Then, we find 
a solution [A'^U)] to give a value A' a (U)-U/2 ~ -1.427 
in the limit. While, due to the lack of effects from the 
higher-order processes in the kinetic energy term, the ap- 
proximate value deviates from the numerical estimation, 
this exhibits the following, i.e., the perturbative expan- 
sion becomes singular on the 2A = U line so that the 
phase boundary deviates from the line. This singularity 
also exists in the perturbative calculations of the SDW 
and CDW state energies (Ssdw and Eqdw)- And then 
the direct transition line between these phases cannot be 
determined from the equation Egnw = -Ecdw, which is 
highly contrasted to the case of the extended Hubbard 









, 


- \au.) ' 

A (U I) 


p , 






+ 8-10, 10 
» 10-12, 12 
a 12-14, 14 
- 14-16, 16 
° 16-18, 18 

1 1 1 


1 1 1 1 


1 1 1 1 





0.5 


1 

u 



FIG. 6: The deviations of boundaries from the 2A = U 
line, A V (U,L) - U/2. We use u = U/(U + 4) as the x axis. 
The correspondence between marks and system sizes is given 
in the figure. Marks with solid (dotted) curves exhibit the 
deviations in the spin (charge) part. 

model (EHM) including the nearest-neighbor Coulomb 
interaction £V Vrijnj + i^2. Since the spin-charge cou- 
pling term with the dimerized spin part generates one of 
the relevant forces, A p (U) should be affected by that of 
the spin part. Besides the present model, it is known that 
EHM possesses the coupling term V cos \fd>4> p cos yl^ 
in its bosonized form^ and that the BCDW state with 
the locking points (Vo^.o-) = is stabilized around the 
2V = U line in the weak- and intermediate-coupling 
region^ The corrections to g v from higher-energy states 
stabilize itj2£ but the coupling term forces the boundaries 
to merge into the single first-order phase transition line 
between the SDW and CDW states in the strong-coupling 
region because it raises the BCDW state energy. How- 
ever, in the present BCDW state, the locking point 4>o m 
Table [I] may take a value so as not to bring about a large 
energy cost due to the coupling term Eq. Therefore, 
the existence of the SDI phase is not prohibited even in 
the strong-coupling limit in contrast to the EHM case. 
Of course, these arguments are qualitative and intuitive 
ones, so an effective theory in this limit is required for 
the precise description on the limiting behaviors. 

Finally, we comment on the Berry phase 
mcthodiSiliiiiiSiiii The Berry phases for the charge 
and the spin parts j u are related to the ground-state ex- 
pectation values of the twist operators as 7^ = Im log z v 
where 

«, = <W>> ^ = (C/tC/; 1 ), (14) 

and U s = exp[(27ri/L) X^=i j n j,s]^ Since z v is real at 
the half filling with zero-magnetic field, 7„ (=0 or n) in- 
dicates the sign of z v . On one hand, z„ can be related 
to the bosonic field as z p , a cx =f(cos v8<^p o-)j so that it 
includes the information of the locking points given in 
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FIG. 7: Behavior of the ground-state expectation value of 
the twist operator z v (y = p, a) near the 2A = U line. The 
correspondence between marks and system sizes are given in 
the figure. 
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FIG. 8: Comparisons of the system-size dependences of the 
transition points obtained by the LC and PRG methods vs 
by the condition z v = 0. The fitting curves show the extrap- 
olations of data to the thermodynamic limit. 



Tabic E 15 In Fig- El we show behaviors of z„ near the 
2A = U line for U = 16 and find that with the increase 
of A both of these increase and change their sign. As 
shown in the lower panel of Fig. |S| the condition z a = 
gives a close value to the result of the LC method, so 
it may provide a proper estimation of the spin-gap tran- 
sition point A CT ii^4i& On the other hand, the zero point 
of z p exhibits a deviation from the PRG result (see the 
upper panel of Fig.[SJ). Since (f>o continuously varies with 
A, z p can take a finite value on the 2D-Ising transition 
point in the thermodynamic limit, which is highly con- 
trasted to z a on the spin-gap transition point. In fact, the 
size-dependent zero points are seemingly extrapolated to 
a value different from our PRG estimation, so that the 
condition z p = does not specify the transition point. 



On the other hand, we also find in Fig. [7| that there is a 
point A ~ 7.3 at which z p is almost independent of L. 
This crossing point is expected to be a good estimator for 
the 2D-Ising transition point in the charge part A p be- 
cause this is quite close to the PRG result even for small 
L. However, a theoretical explanation of this possibility 
is still open. 

To summarize, we have investigated the ground-state 
phase diagram of the one-dimensional half-filled Hubbard 
model with the alternating potential, especially in order 
to verify the scenario given by Fabrizio, Gogolin, and 
Nersesyan, we have numerically treated the phase tran- 
sitions observed in the spin and charge parts: We cal- 
culated the spin-gap transition points A CT in the spin 
part by the level-crossing method (see also the argu- 
ment for the spin-gap transition in Ref. Il6|) and the two- 
dimensional Ising transition points A p in the charge part 
by the phcnomenological rcnormalization-group method. 
We confirmed that, adding to the Mott and band insula- 
tors, the "spontaneously dimerized insulator" accompa- 
nied by the long-range-ordered 2fcp bond charge-density 
wave is stabilized as the intermediate phase for all U > 0. 
Then we checked the SU(2)-symmetric Gaussian (2D- 
Ising) criticality of the spin (charge) part by treating the 
low-lying excitation levels in the finite-size systems, and, 
simultaneously, wc performed the finite-size-scaling anal- 
ysis of the charge-excitation gap to clarify the critical 
phenomena around A p . The comparison with the relat- 
ing work was performed to check the reliability of our 
numerical results and to exhibit the efficiency of our ap- 
proach. 

After submission of this paper, we became aware of the 
work investigating the ground-state phase diagram and 
the universality of the transition in the charge part by 
the use of finite-size-scaling analysis of the DMRG calcu- 
lation datarfS They have found two transition points and 
succeeded to obtain v = 1 in agreement with our conclu- 
sion, while the estimated exponent for the susceptibility 
of the BCDW order parameter shows a deviation from 
the theoretical value 771 = 1/4, e.g., r\\ ~ 0.45 at the point 
on the BI-SDI phase boundary A = 10 and U cl = 21.385 
(in their notation) . In this paper wc have treated the el- 
ementary excitations in the TLL system specified by the 
discrete symmetries of the lattice Hamiltonian with the 
twisted boundary condition, whereas they have measured 
the BCDW order parameter, (i.e., a composite excitation 
of the spin and charge degrees of freedom) with the larger 
energy scale. 
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